Dual Fermion Dynamical Cluster Approach for Strongly Correlated Systems 
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We have designed a new multi-scale approach for strongly correlated systems by combining the 
Dynamical Cluster Approximation (DCA) and the recently introduced dual-fermion formalism. This 
approach employs an exact mapping from a real lattice to a DCA cluster of linear size L c embedded 
in a dual fermion lattice. Short-length-scale physics is addressed by the DCA cluster calculation, 
while longer-length-scale physics is addressed diagrammatically using dual fermions. The bare and 
dressed dual fermionic Green functions scale as 0(1/L C ), so perturbation theory on the dual lattice 
converges very quickly. E.g., the dual Fermion self-energy calculated with simple second order 
perturbation theory is of order 0(1/ L^), with third order and three body corrections down by an 
additional factor of 0(1/L C ). 



I. INTRODUCTION 

Dynamical mean-field theory-ir— has been remarkably 
successful at capturing the physics of strongly correlated 
systems dominated by spatially local correlations. Suc- 
cesses include the description of the Mott transition in 
the Hubbard model, screening effects in the periodic An- 
derson model, as well as the description of correlation 
effects in realistic systems^— . 

Since the introduction of the Dynamical Mean-Field 
Approximation (DMFA) there have been a number of at- 
tempts to develop formal extensions around the DMFA 
that incorporate non-local corrections. These include 
cluster extensions of the DMFA, such as the Dynamical 
Cluster Approximation (DCA)2r— or the Cellular Dynam- 
ical Mean-Field Theory (CDMFT)^, and multi-scale ap- 
proximations where the DMFA or DCA vertices are used 
to parameterize two-particle field theories and longer 
ranged correlations can thus be captured^— . One of 
the main limitations of these theories is that they con- 
verge slowly with the linear cluster size L c , especially for 
the calculation of transition temperatures. 

The Dual Fermion formalism^— is however, dis- 
tinctly different from other cluster extensions of the 
DMFA. In the dual fermion formalism, the lattice action 
is first mapped onto a dual fermion action where the in- 
teraction vertices are the n-body reducible vertices of the 
cluster. This mapping is exact, so the dual fermion for- 
malism provides a complete and exact formalism for the 
lattice problem. Thus far, the dual fermion formalism 
has only been explored using the DMFA or the CDMFT 
as cluster solvers^. However, the CDMFT has the dis- 
advantage in this context that it violates translational 
invariance, so that the CDMFT vertices are rank-4 ten- 
sors in the spatial or momentum indices, which are too 
large to be stored and manipulated on most computers, 
especially for large clusters. Thus in this manuscript we 
propose the Dual Fermion Dynamical Cluster approach 
(DFDCA) , within which the long-ranged correlations can 
be systematically incorporated through the dual fermion 
lattice calculation. Since the DCA preserves the trans- 



lational invariance of the lattice system, the DCA two- 
body vertices are rank-3 tensors which, for modest clus- 
ter sizes, will fit in the memory of modern computers. 
Another difference, which we will discuss in detail, is 
that the small parameter for the DFDCA is the dual 
fermion single-particle Green function, which scales as 
Gd ~ 0(1/ L c ) with L c being the linear cluster size. As 
a result, perturbation theory on the dual fermion lattice 
converges very quickly. Simple second order perturba- 
tion theory on the dual fermion lattice already yields a 
dual fermion self-energy of order 0(1/ L^/) with two-body 
and three-body corrections down by an additional factor 
of 0(1/L C ). Higher order approximations are also possi- 
ble, since the, e.g., three-body vertex corrections to the 
DFDCA self-energy are small, 0(1/ Lf). Therefore, the 
resulting DFDCA formalism converges very quickly with 
increasing cluster size, with corrections to the self-energy 
no larger than 0(1/L^). 



II. MAPPING THE DCA FORMALISM TO 
DUAL FERMIONS 

We will derive the DFDCA formalism with the exam- 
ple of the Hubbard model. Its Hamiltonian is 



H = - *«j(4r c 3* +h.c), 



<ij> 



C/^(n it -l/2)K;-l/2) 



(1) 



where tij is the matrix of hopping integrals, c^J is the 
annihilation (creation) operator for electrons on lattice 
site i with spin a, xii a = c^ ia Ci a , and U the intra-atomic 
repulsion. 

The DMFA, and its cluster extensions such as the 
DCA, are based upon the common idea of embedding 
a cluster in a lattice. We assume that the cluster, of 
size N c = L®, dimensionality D, sites labeled by I and 
wavevectors K, is embedded in a large but finite-sized 
lattice of size N with sites i and wavevectors k. In the 
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DCA, the reciprocal space of the lattice is divided into 
iV c cells of identical geometry and linear size Afc. The 
cell centers are labeled by K, and the points surrounding 
K within the coarse-graining cell are labeled with fc. We 
will also invoke a dual space lattice which is of the same 
size and geometry as the real lattice. 
The action for this model is 

S[c*,c} = - ^2 c l,k,A( ioj +l- L ) 1 - h k] c ^,k, a +^2Si oc [c*,Ci}, 

(2) 

where Si oc [c*,Ci] is the local part of the action including 
the Hubbard interaction term, c* and Ci are now Grass- 
mann numbers corresponding to creation and annihila- 
tion operators on the lattice, fj, the chemical potential, 
hk the lattice bare dispersion, and u> = (2n + 1)ttT the 
Matsubara frequencies. Decomposing the wavevector ac- 
cording to k = K + k, the lattice action becomes 

S[C*,C] = ]P Sl oc [c* , Ci] 

i 

- J2 c t,K+kJ^ + l^ 1 ~ h K+kK,K+k,A 2 ') 

uj,K,k,a 

The goal is to express this action in terms of the DCA 
cluster problem^ 

•^cluster I s * ' c l = /,Sloc&,ci] 
i 

where c*j and c/ are now Grassmann numbers correspond- 
ing to creation and annihilation operators on the DCA 
cluster, and A(-K", iui) is the cluster hybridization func- 
tion. To this end, we add and subtract the hybridization 
function and coarse-grained dispersion, i.e., 

Ld,K,k,(T 

= Jf 2^ C*,,K,A hK + A ( K , iul )]c Wt K t „, (5) 

where the last line follows from the DCA coarse-graining 
identity 

Ki,K,^,K,a = -Jj ^2 C u,K+k,<T C u,K+k,<r ( 6 ) 
k 

and the coarse-grained dispersion is given by 

^ K = ^-H h K+k- ( 7 ) 

fc 

The DCA coarse-graining identity preserves the 
Fcrmionic Lie algebra, despite the fact that it is not a 
canonical transformation, 

fc 

(8) 



where the last step follows since the coarse graining cells 
surrounding K and K' have the same geometry and con- 
tain the same number of states which, therefore, may be 
labeled with the same fc. We then obtain 

S[c*,c] = ^2Si oc [c*,Ci] 

i 

~ ^.K+kJ^ + ^-^K- A{K,iu>)}c uK+ ~ ka 

Ld,K,k,(7 

- c t,k,AMM(k),iu) + h M{k) - h k ]c^ K<T . (9) 

uj,k,a 

In the third line of this equation we have introduced the 
function M (fc) which maps the momentum fc in the DCA 
momentum cell to the cluster momentum contained in 
that cell. Coarse-graining the first and the second terms 
on the right hand side of the above equation yields the 
cluster action Since the latter is independent of fc, 
we may write 

S[c*,c) = J2 ^cluster t 5 *> 5 ] 
fc 

- £ c*, fe , CT [A(M(fc),iu;) + ^M(fc) — h k ]c u .k,a- (10) 

D,k,a 

Again, up to this point, we have only re-arranged terms 
and employed an identity which defines c. No approxi- 
mation has been made. 

The dual fcrmions are now introduced by means of the 
following Gaussian identity 

J cx P(-/>«/3//3 ~ flbapcp - dZ t b a /}fp)ILydf*cff- l 

= det(a)exp[c* {ba~ 1 b) a pcp\ (11) 

for Grassmann variables in the path integral representa- 
tion for the partition function 

j eM-S[c*,c})V[c*,c]. (12) 

To be specific, we choose the (diagonal) matrices accord- 
ing to 

acj,k,a = g~ 2 {M(k),iu>)[A(M(k),iui) + h M(k) - h^ 1 ; 
b u ,k,a = rHM(k),iu). (13) 

where g is the single particle Green function calculated 
on the DCA cluster. Applying the above identity to the 
second term in (fTU|) yields 

J w,k,a fu,k,<r 

fa o g 2 (M(k), iu)[A(M(k), iu>) + h M {k) - h k ] 
+ E \fZ,k,*r\M{k)M<>»,k* + h.c]. (14) 

aj,fc,cr 

The essential observation now is that, since 
<7 (ilf(fc), iui) = g~ x (K, iu>) is independent of fc, 



3 



the second line of (fT4]l may be coarse-grained using 
again the DCA coarse-graining identity 



(15) 



As a consequence the lattice action, Eq. (fTO)) . can be 
expressed in the form 



S[c*, c; /*,/] = J2 S TCstT [c*,c-j\.f} 

+ E 



f* f 

J u:,K+k,cr Ju,K+k,o 



g 2 {K, iu)[A(K, iu) +h K - h k ] 



(16) 



where 

^rcstr [c* > c; / , /] = Cluster [ s * > 5 1 

is the action restricted to the cluster. 

The transformation to dual fermions is completed by 
integrating out the fermionic degrees of freedom corre- 
sponding to c and c*. Since S^estr l& independent of k, 
this can be done individually for each cluster 



1 



exp(-S restr [c*,c;/ ,f])V[c*,c] 



cluster 

= cx p f - E T u ,K,*r\KMJ»,K, a - v[TJ] 



,K,i 



(18) 



Eq. (|18|1 defines the dual potential which can be obtained 
by expanding both sides and comparing the resulting ex- 
pressions order by order. It is given byi^: 

vF,7\ = \ E E E 

KK'Q UJU)' V (71 ,(T2,(T3,CT4 

7(T U (T2,<r3,(T4.(K, K', Q; iu, iu' ', iv) 

x f ui+v,K+Q,ai fui,K,a2 fu)',K',<r a fu'+v,K'+Q,Ui 

+ ... (19) 

where 7 is the full (reducible) vertex of the cluster quan- 
tum impurity model, and the higher order contributions 
involve the n-body (for n > 2) reducible vertices as the 
bare interaction. Integrating out the lattice fermions re- 
sults in an action which depends only on the dual fcrmion 
degrees of freedom given by 

Sd[f*,f] = -Yf: ka G° d (k,iu)- 1 u„ + J2v[T,7l?o) 

where G® d is the bare dual Green function defined by 

g(K,iu) 2 



G° d (k,iu) = - 



g(K, iu) + (A (if, iu) + h K - h k ) ' 



This quantity together with the dual potential V[f , /] 
provides sufficient input for a many-body diagrammatic 
perturbation calculation on the dual lattice. 

Note that besides the DCA coarse-graining process in- 
troduced here, the above derivation is a natural general- 
ization of the dual fcrmion DMFA formulation of Rubtsov 
etalM 



A. Self-consistency condition 

In rewriting the lattice action in terms of the cluster 
impurity model in the above derivation, the DCA hy- 
bridization function has been added and subtracted and 
hence is an arbitrary quantity. In order to fix this quan- 
tity we impose the condition 

G° d (K, iu) = ^ £ G d (K + k, iu) I 0. (22) 

k 

To appreciate the consequences of this condition, first 
consider the DCA lattice Green function 

G^ CA (K+k,iu) = (iu+n)l-h K+k —'Z c (K,iu), (23) 

which can be expressed in terms of the cluster Green 
function 

g (K, iu) = (iu + fj,)l — hfc — T* C (K, iu) — A C (K, iu), 

(24) 

as 



Gj^ CA (K+k,iu)=g 1 (K,iu) + A c (K,iuj) + h K -h K+ ~ k . 

(25) 

Using the last expression, one may straightforwardly de- 
rive the following identity relating the DCA lattice Green 
function to the bare dual Green function 

G d '°{K + k,iu) = G BCA (K + k,iu) -g(K,iu). (26) 



(21) 



Hence the above condition (|22|) is equivalent to requir- 
ing that the coarse-grained DCA lattice Green function 
be equal to the Green function of the cluster impurity 
model. This is exactly the DCA self-consistency condi- 
tion. The DCA solution is therefore obtained if no di- 
agrammatic corrections are taken into account and the 
hybridization is determined such that (|22|) holds. Conse- 
quently, we have a perturbation theory around the DCA 
as the starting point. While the DCA solution only de- 
pends on the cluster momentum K, the dependence on k 
can be introduced by solving the dual problem pcrturba- 
tivcly. 



B. Scaling of the Dual Fermion DCA approach 
with cluster size 



The bare dual Green function is given by Eq. (|2ip . 
If we introduce the linear cluster size L c through N c = 
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L®, one finds that the term (A(K, iu) + ha — hk) ~ 
0(1/ L c ). The small nature of this term for large L c 
should ensure rapid convergence of the DFDCA. In par- 
ticular, we then have 

G° d (k, iu) = -g(K, iu) (A(K, iu) + h K - h k ) g(K, iu) 
+0(1/L% (27) 

i.e., the bare dual Green function also scales like 

G d (k,iu)~O(l/L c ). (28) 

However, at points of high symmetry, where hx — h k ~ 
0(1/L 2 C ), G° d (k,iu) will fall more quickly than 0(1/L C ). 

To illustrate the typical scaling behavior of the bare 
dual Green function, we plot in Fig. [TJ as a function of 
1/L C , the ratio of \G d (k,iu = iirT,L c )\ averaged over k 
to the average \G d (k,iu = iirT,L c = 1)|. We also plot 
the average of the ratios. The former initially falls more 
quickly than 0(1/ L c ), while the latter displays a slower 
initial slope. However, for L c > 4 both fall roughly lin- 
early in 1/L C . This behavior is found to be independent 
of temperature (not shown) , since it is a purely algebraic 
effect. 




FIG. 1: Scaling plot for the bare dual Green function. Here 
we have used the 1-D Hubbard model to analyze its scaling 
behavior. Except for very small L c values, the two ratios scale 
linearly according to Eq. (1 281) . 

Applying the standard tools to the dual fcrmion action, 
one obtains the formal expression 

G d (k, iu) = G° d (k, iu) + G° d (k, iu)T d (k, iu)G° d (k, iu), 

(29) 

for the full dual fcrmion Green function G d (k, iu), where 
the reducible self-energy or scattering matrix T d (k, iu) of 
the dual system is introduced. We will show later, that 
T d (k,iu) will be at most of order 0(1/ Lfy, and we can 
infer the scaling 

G d (k,iu)~0(l/L c ) (30) 
for the full dual fcrmion Green function too. 



Once the dual fcrmion Green function is known, one 
can reconstruct the real lattice Green function as 

G(k,iu) = g(K,iu)~ 2 (A(K,iuj) + h K - h k )~ 2 G d (k,iu) 

+ {A(K,iu) + h K -h k )~ 1 . (31) 

Since G(k, iu) is the Green function of the real lattice, it 
should scale as 

G(k, iu) ~ 0(1) (32) 

with respect to any length scale. On the other hand, for 
the two terms on the right hand side in (|3"Tj) we find 

g(K, iu)- 2 (A(K, iu) +h K - h k ) ~ 2 G d (k, iu) ~ 0(L C ) 

(33) 

and 

(A(K, iu) +h K -h k )~ 1 ~ 0(L C ) . (34) 

Thus, the two 0(L C ) terms must cancel each other. To 
verify this requirement, we insert the zeroth order contri- 
bution of the dual Green function into the original Green 
function, and after some algebra we indeed obtain 

G(k, iu) ~ g(K, iu) - 0(1), (35) 

with a correction given by 

AG(k,iu) ~T d (k,iu). (36) 

Therefore, the correction to the real Green function 
through the dual fermion approach scales the same way 
as the dual self-energy. 

Presently, the dual potential Eq. (flH|) still contains an 
infinite hierarchy of vertices. The previous discussion 
now provides a very important insight into the contribu- 
tions of these vertices to a perturbative expansion: Each 
?i-body diagrammatic insertion will involve a vertex and 
n Green function lines. In the parameter region away 
from a critical point the dual potential will be of order 
0(1). As noted before, the dual Green function is of or- 
der 0(1/L C ), i.e., each dual space diagrammatic insertion 
is of order 0(l/i^) when it involves the two-body dual 
space interaction, of order 0(1/Lj?) for the three-body 
interaction, and so on. This means that the two-body 
contribution to V, explicitly shown in Eq. (fT^)l . will ac- 
tually dominate and low-order perturbation theory will 
be sufficient to accurately capture the corrections to the 
DCA from the dual fermion lattice. 



C. Mapping back from the Dual-Fermion to the 
Real Lattice 

The relation of the real fermion Green function to the 
dual Green function has been been established in Eq. 
([3Tj) . This is an exact relation which follows by taking 
the functional derivative of two equivalent partition func- 
tions. They are linked through the same Gaussian iden- 
tity that has been used to introduce the dual fcrmions 
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(Eq. (jTTjl ). Higher order derivatives then allow us to 
derive relations between higher order cumulants. From 
this recipe, we find the following relation between the 
two-particle reducible vertex functions 

F k ,k',q;iu,iu;',i V = T(k + q,iu + iv)T(k,iuj) F^ k , tq;iu , iu > >iv 
x T(k' ,iu)')T(k' + q,iu' + iv) (37) 

in real and dual space, where 



G(k, iu)(A(K, iuj) + h K - h k )g(K, iui) 
= -[l + g(K,iw)H d (k,iu))\~ X - (38) 

Similar relations hold for many-particle vertex functions. 
With the help of the two-particle vertex function we 
can now express the corresponding susceptibility as x = 
Xo + Xo^Xo- Since from Eq. ([38} it follows that T(k, iuj) 
is always finite, a divergence of x, signaling an instability 
or phase transition in real space, necessarily corresponds 
to an instability in the quantity F d in the dual fermion 
space. In order to locate the instabilities, it is hence suf- 
ficient to search for a divergence of the Bethe-Salpeter 
equation in the dual space. For the special case when 
no diagrammatic corrections to the dual self-energy and 
vertex are taken into account, T(k,iu) = —1 and both 
DFDCA and DCA would produce the same phase di- 
agram. In general cases, the DFDCA will produce re- 
sults more realistic than DCA due to the inclusion of ad- 
ditional long-ranged correlations from the dual fermion 
lattice diagrammatic calculation. 

III. DUAL FERMION DIAGRAMS 

In the DFDCA formalism, the dual fermion Green 
function is 0(1/L c )(c.f. Eq. ([28]) ). i.e., it acts as the small 
parameter in the diagrammatic expressions. In addition, 
in the strong coupling limit, the Green function is pro- 
portional to the hopping t/U r^- so each Green function 
leg contributes a factor of 0((t/U) / L c ). In the weak cou- 
pling limit, the Green function remains C(l/L c ), but the 
vertices are now small, with the two-body vertex behav- 
ing like 0(U/i), the three-body vertex like 0(U 2 ), and 
so on. Each two-body diagrammatic insertion, composed 
of a two-body vertex and two dual fermion Green func- 
tion legs, then scales like 0(1/L 2 ), with an additional 
factor of U or t 2 in the weak and strong coupling limits, 
respectively. Each three-body diagrammatic insertion, 
composed of a three-body vertex and three dual fermion 
Green function legs, scales like 0(1/ L 3 ) with additional 
factors U 2 or t 3 in the weak and strong coupling limits, 
respectively. 

The boundary condition Eq. ([2"2"j) also constraints the 
diagrammatics. For example, the first-order contribution 
to the dual self-energy from the 2-body interaction is the 
Hartree-Fock contribution shown in Fig. [5J Since the 
vertex depends only upon the small cluster K, the dual 




FIG. 2: Lowest order contributions to the dual fermion self- 
energy from the two-body interaction (left) and to the two- 
body interaction from the three-body term (right). Since the 
bare n-body vertices depend only upon the small cluster K, 
the dual Green function line may be coarse-grained and are 
therefore zero according to Eq. (|22l) . 

Green function line may be coarse-grained. The result 
is zero by virtue of Eq. (|22|) . Physically, this term must 
be zero since the Hartree term is already included in the 
cluster contribution to the self-energy. Therefore, the 
first finite contribution to the dual fermion self-energy 
comes from the second order graph which contains three 
dual-fermion Green function lines. This and all higher 
order contributions described by the Schwinger-Dyson 
equation already arc of order 0(1/L 3 ), or smaller. There- 
fore, the fully dressed dual fermion Green function retains 
the scaling of the bare dual fermion Green function 

G d (k,ioj) ~C7(1/L C ) (39) 

as anticipated earlier. 

Similarly, the first-order 3-body contribution to the 
dual two-body vertex, also shown in Fig. [21 is zero. To 
see this, suppose the top leg is labeled by momentum 
k = K + k. Since the remainder of the 3-body vertex 
does not depend upon k, we may freely sum over this 
label. Again, the result is then zero through Eq. (|22|). 

As the cluster size becomes large, the DFDCA clus- 
ter problem may be accurately solved using low order 
perturbation theory, keeping only the 2-body interaction 
vertex. As described above, two-body vertex insertion 
contributes an extra factor of 0(1/L 2 ), while three-body 
vertex insertion contributes an extra factor of 0(1/ L^). 
It is therefore possible to use standard perturbation the- 
ory based on a two-body vertex to solve the dual-fermion 
DCA cluster problem, with an accuracy which turns out 
to be at least of order 0(1/L^). 

For example, simple second order perturbation the- 
ory yields a self-energy 0(1/ L^). Two-body corrections, 
composed of a two-body vertex and two further Green 
function legs will contribute an extra factor 0(1/ L 2 ) . 
The first three-body contribution is the second order 
graph composed of one 2-body vertex and one 3-body 
vertex. It has four internal Green function legs, and 
is of order 0(1/Lf) so that the first three-body correc- 
tion is smaller than the simple second order dual fermion 
self-energy composed of 2-body vertices by a factor of 
0(1/ L c ). Self consistency, needed to impose the bound- 
ary condition Eq. (|22[) . is more important for the self- 
energy than higher order or three-body contributions. 

Generally, the leading non-trival n-body (n > 3) ver- 
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tcx contribution to the self-energy is constructed from 
an n-body vertex and an (n — l)-body vertex, which are 
connected by (2n — 2) internal legs, as shown in Fig. [3] 
Thus this contribution scales as 0(1/ L%. n ~ 2 ). 



$ = $ 




2n-2^ 



FIG. 3: Leading non-trivial n-body vertex contribution to 
the self-energy. It is contracted by one n-body vertex and 
one (n-l)-body vertex. Since there are (2n-2) internal legs, 
this contribution scales as 0(1/L^. n ~ 2 ). 

As another example, consider the equation for a tran- 
sition, in the pairing matrix formalism (Fig. [4]) 



r dX > = $ 



(40) 



where $ is the leading eigenvector of the pairing matrix 
TdX d - A transition is indicated by the corresponding 
eigenvalue approaching one. To lowest order, the irre- 
ducible dual fermion vertex « 7 is just the bare dual 
fcrmion interaction, and the legs in Xd are n °t dressed 
by the dual fcrmion self-energy. In this case the transi- 
tion temperatures of the DCA are reproduced (e.g., sec 
Fig. [SJ. The lowest order corrections to the DCA come 
from the second order corrections to the vertex, which 
contain two dual fermion Green function legs and are 
therefore 0(1/ L 2 ). The low order contributions to Xd 
are x° d ~ G di l + ^dG° d + ■ ■ -)G d , and thus the lowest 
relative correction to Xd, ^ s °f or der 0(1/ L 4 ). Therefore, 
the cross channel second order corrections to the vertex 
are more important than the second-order corrections to 
the self-energy when the DCA cluster size is large. We 
note that this is not only true for the DFDCA, but also 
for the DFDMFA in the strong coupling limit where the 
small parameter t/U replaces l/L c . Furthermore, higher 
order approximations such as the ladder approximation 
that do not include these cross channel contributions arc 
not appropriate for the solution of the dual fcrmion lat- 
tice in the limit of large DCA cluster size or small t/U. 

Higher order approximations like the fluctuation- 
exchange approximation (FLEX)2£, which include the 
cross channel contributions to T d , should on the other 
hand be quite accurate. In fact, the FLEX contains all 
two-particle diagrams to third order. The first diagram 
neglected by the FLEX is composed of one three-body 




owl) 




owi) 



FIG. 4: (a) Equation for T c . Transition temperatures on the 
dual fermion lattice are identical to those calculated on the 
real lattice, (b) The low order corrections to the dual fermion 
irreducible vertex Fd- The second order terms are of the or- 
der 0(1/1%) with corrections 0(1/1%). (c) Contributions to 
the dual fermion self-energy. It is dominated by the second- 
order term which is of the order 0(1/ L^.) with corrections 
0(1/ L^). The self-energy adds relative corrections to Xd °f 
order 0(1/ L^) (see the text for the detail), so the most im- 
portant contributions to the equation for T c come from the 
second-order cross channel contributions to IV 



and one two-body vertex and would contribute a correc- 
tion 0(1/ L 4 C ) to the self-energy or 0(1/1%) to the ver- 
tices. 



IV. RESULTS 

In this section, we will present numerical results from 
a DFDCA calculation, where the interaction expansion 
continous-time quantum Monte Carlo method^ is em- 
ployed to solve the cluster problem within the DCA cal- 
culation. We will restrict the discussion to the two- 
dimensional Hubbard model on the square lattice with 
only nearest neighbor hopping. Thus, for half- filling 
we expect strong antiferromagnetic correlations, which 
will drive an antiferromagnetic transition within DCA. 
In this the Mermin- Wagner theorem prohibits 

long-range order except at zero temperature, we expect 
strong renormalization of the Neel temperature, T/v, from 
DFDCA. 

To check the correctness of our implementation of the 
DFDCA approach, we first carry out calculations with 
the correction from the dual-fermion lattice turned off. 
For this trivial case, one expects DFDCA to reproduce 
the same physics as DCA. Fig. [5] displays the leading 
eigenvalues for different cluster sizes at filling (n) =0.95 
for the antiferromagnetic channel. Note that for each 
cluster size, both the DFDCA and the DCA leading 
eigenvalues cross the line A = 1 at the same tempera- 
ture, which is the mean-field Neel temperature, and that 
with increasing cluster size Tn decreases, as expected. 
It is also interesting to note that the DFDCA provides 
a sensitive way to monitor the finite-temperature transi- 
tions since the DFDCA leading eigenvalues have a steeper 
slope when crossing the A = I line. 

For the non-trivial DFDCA calculation, we expect to 
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T/4t 



FIG. 5: (Color online) Plots of leading eigenvalues for dif- 
ferent cluster sizes for the anti-ferromagnetic channel with 
U = 6t and filling (n) — 0.95. Lines without symbols are re- 
sults from the DCA calculation for clusters with sizes N c = 8, 
12 and 16, while lines with symbols are results from the 
DFDCA calculation without self-energy correction. For the 
latter, we have used a linear size of the dual fermion lattice 
as large as several hundreds (JV = L x L, L ~ 200). The inset 
is an enlarged view around the transition point. Note that 
both calculations produce the same transition temperatures 
as expected. 
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FIG. 6: (Color online) Plots of leading eigenvalues for the 
anti-ferromagnetic channel. They are calculated with differ- 
ent approximate methods in the dual-fermion lattice. The 
parameters used are U = 4t, (n) — 1, the DCA cluster of size 
N c = 1, and the dual fermion lattice of a size iV = 4 x 4. 



see for a fixed cluster size a reduction of the Neel tem- 
perature since correlations beyond the cluster scale are 
now incorporated by the dual fermion calculation. For 
the dual fermion lattice, we employ different approxima- 
tion schemes: the self-consistent second-order perturba- 
tion theory (SOPT), FLEX and the parquet approxima- 
tion (PA)2i. The results are collected in Fig. [51 where 
the power of DFDCA manifests itself clearly. The sim- 
ple second-order correction from the self-energy is already 
able to reduce the Neel temperature by ten percent. Tak- 
ing into account more Feynman diagrams with higher or- 
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FIG. 7: (Color online) Plots of the inverse anti-ferromagnetic 
and d-wave pairing susceptibilities calculated with different 
approximate methods in the dual fermion calculation. The 
parameters used are U = 8t and N c — 4. The linear depen- 
dence of the results with 1/L 2 (see Fig. [8} is used to extrap- 
olate the L — oo limit results. 



ders, for example by FLEX or PA, continues to reduce 
the Neel temperature. However, the inclusion of vertex 
correction tends to increase the Neel temperature again. 
For example, the eigenvalues labeled Eflex are calcu- 
lated with a bare dual fermion vertex and FLEX dressed 
legs, while those labeled Eflex + Tflex are calculated 
with both FLEX dressed legs and vertex (see b and c in 
Fig. [4] for contributions up to second-order in the bare 
dual fermion vertex 7). 

Up to now we have only discussed the leading eigenval- 
ues of the vertex function. Of course, the DFDCA also 
allows to calculate the full susceptibility from the Bethe- 
Salpetcr equation. Two typical results for an N c = 2 x 2 
DCA cluster are shown in Fig. [3 as function of temper- 
ature for U = 8t. In the left panel, the inverse stag- 
gered susceptibility for half filling is displayed, while the 
right one contains results for the inverse d-wave pairing 
susceptibility at a filling (n) = 0.95. Due to the heavy 
computational cost for the parquet calculation, here we 
only use the SOPT and FLEX in our dual fermion lat- 
tice calculation. Although L c = 2 is not really large, the 
DFDCA is still able to significantly reduce the mean-field 
Neel and the abnormally large superconducting transi- 
tion temperatures. It is quite interesting to note, that 
for the anti-ferromagnetic channel at half-filling, SOPT 
and FLEX produce similar results, both being different 
from the DCA results. The effect of vertex correction is 
small in this case. For the d-wave pairing susceptibility 
at (n) = 0.95, on the other hand, SOPT in dual space 
makes almost no difference from the DCA results, but 
the FLEX tends to significantly reduce the pairing sus- 
ceptibility. Again, the inclusion of vertex correction has 
the opposite effect, i.e., it leads to a slight increase of the 
critical temperatures. 

In the derivation of the DFDCA approach, we have 
assumed that the dual fermion lattice is infinite. How- 
ever, in practical calculations, the size is limited due to 
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FIG. 8: (Color online) The L dependence of the leading eigen- 
value for different DCA clusters. The parameters used are 
U = 8t, (n) = 1, T, dual = and /3 = At/T. The nice lin- 
ear dependence of 1/L 2 can be readily observed, which is due 
to the periodic boundary conditions used in the dual-fermion 
calculation. 



the algebraic increase of the computational cost. This 
results in some deviations from the infinite size system. 
Fig. [5] shows the L (N = L x L) dependence of the lead- 
ing eigenvalues for different DCA clusters. The nice lin- 
ear dependence of the leading eigenvalues on 1/L 2 can 
be readily observed. This is due to the periodic bound- 
ary conditions used in the dual-fermion calculation. This 
property allows us to reduce the computational cost of 
our calculation by using two small L's and extrapolating 
to obtain a rather accurate approximation of the L = oo 
result. 



DISCUSSION 



QMC simulations for larger clusters, and access a wider 
region of parameter space. 

The DFDCA has an additional advantage that it 
is parameterized by the full (reducible) vertex func- 
tion calculated on the DCA cluster. Other multi-scale 
method*^— rely upon the calculation of the cluster ir- 
reducible or fully irreducible vertices. Our recent numer- 
ical experiments show that inverting the Bcthc-Salpctcr 
equation to obtain the irreducible vertex, which is also 
the first step in the calculation of the fully irreducible ver- 
tex, fails in some parameter regions, especially for large 
U or near half filling. This difficulty is avoided in the 
DFDCA. 

The dual fermion mapping also assumes that the dual 
fcrmions are treated on an infinitely large lattice. In 
practice however, they are treated on a finite-size lat- 
tice. Thanks to the finite-size scaling behavior observed 
in Fig. [51 finite-size calculations are used to extrapolate 
to the infinite-size lattice, leading to a reduction of the 
computational cost in the dual fermion lattice calcula- 
tion. 

Note that in the calculations presented here, we have 
not performed the full self-consistency where the dual 
fermion result is used to determine the DCA cluster hy- 
bridization that is fed back into the DCA calculation un- 
til convergence. However, this first iteration already pro- 
duces more satisfactory values for the Neel temperature 
as well as the d-wave superconducting transition tem- 
perature. We can anticipate that the full self-consistency 
will further improve the performance of this approach. 

With the full self-consistency implemented and tested, 
we are planning to apply this approach to map out the 
phase diagram for the 2-D Hubbard model in the hole- 
doped region. In addition, we are also working on apply- 
ing this approach to the Falicov-Kimball model and the 
Anderson disorder model recently. 



The dual fermion mapping as discussed in section |TT] 
is exact, and the approximation is made only when 
performing the diagrammatic calculation for the dual 
fermion lattice. Justified by the scaling behavior of the 
dual fermion Green function, it suffices to consider the 
2-body term of the interaction and use low order per- 
turbation theory. Correlations beyond the DCA cluster 
size are systematically restored through the dual fermion 
calculation on the lattice. In this sense, the DFDCA 
can be seen as a diagrammatic expansion around DCA. 
This is manifested in Fig. [5] and [6] where we see that not 
including the dual fermion self-energy and vertex correc- 
tions reproduces the DCA transition temperature. When 
these corrections are included, we observe a systematic 
suppression of the DCA transition temperature resulting 
in a more realistic value. This is clearly seen in Fig. [6] 
and [JJ Since correlations at intermediate length-scale 
are taken into account by the dual fermion lattice calcu- 
lation, we can use small clusters in the underlying DCA 
calculation. As a result, we arc able to greatly reduce the 
adverse effect of the minus sign problem encountered in 



VI. CONCLUSION 

We have designed a new multi-scale many body ap- 
proach, the Dual Fermion Dynamical Cluster approach 
(DFDCA), by combining the DCA and the recently intro- 
duced dual fermion formalism. The DFDCA uses both 
single and two particle quantities calculated in DCA as 
the input for the dual fermion calculation. Different self- 
consistent diagrammatic approximations can be used in 
the dual fermion lattice, which systematically restores 
the long-ranged correlation ignored during the DCA cal- 
culation. 

This approach is a systematic expansion around the 
DCA calculation. Our numerical experiments show that 
the zeroth order result (E^ = 0) reproduces the origi- 
nal DCA, and for any non-trivial dual fermion calcula- 
tion, it is an improvement on the DCA calculation. We 
applied different self-consistent diagrammatic methods, 
self-consistent 2nd-ordcr, FLEX and parquet approxima- 
tion, on the dual fermion lattice. They all improved the 
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DCA calculation by reducing the mean-field Neel temper- 
ature by different amounts. In addition, the abnormally 
large superconducting transition temperature of the four 
site cluster calculation can be reduced by this approach 
as well. 
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